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A physical implementation of the adiabatic Grover search is theoretically investigated in a system 
of N identical three-level atoms trapped in a single mode cavity. Some of the atoms are marked 
through the presence of an energy gap between their two ground states. The search is controlled by 
two partially delayed lasers which allow a deterministic adiabatic transfer from an initially entangled 
state to the marked states. Pulse schemes are proposed to satisfy the Grover speedup either exactly 
or approximately, and the success rate of the search is calculated. 
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I. INTRODUCTION 



Search problems can be expressed as finding a set of 
marked items in an unsorted list. The Grover algorithm 
[l| achieves this task quadratically faster than a classi- 
cal algorithm which examines items one by one. It has 
become a paradigm of quantum computation based on 
quantum circuits. The case of two qubits, corresponding 
to a four-element search, has been demonstrated experi- 
mentally in various settings. NMR experiments [2, lla, |4| 
used the spin states of 1 H and/or 13 C in a magnetic field 
as qubits while radio-frequency fields and spin-spin cou- 
plings between the nuclei were used to implement the 
quantum logic gates. In linear optics techniques, the 
individual qubits are represented by different polariza- 
tion or spatial-mode degrees of freedom while the com- 
putation is achieved through essentially a complicated 
interferometer Q or a one-way quantum computer Q. 
Trapped ions systems Q, which present the advantage 
of being scalable, rely on a number of optical and mi- 
crowaves sources to control, entangle and measure the 
qubits which are represented by the ground state hyper- 
fine levels of two trapped atomic ions. There have also 
been proposals of experimental implementations using 
cavity QED where the quantum gate dynamics is pro- 
vided by a cavity-assisted collision or by a strong res- 
onant classical field @. 

A different approach to quantum computation con- 
sists in the controlled evolution of a system obeying the 
Schrodinger equation with a Hamiltonian designed to 
solve a specified problem. It was pioneered by Farhi and 
Gutmann [1 01 ] who considered a time continuous version 
of the Grover algorithm. A given Hamiltonian features 
a marked state whose energy differs from that of the 
N — 1 unmarked ones. A driving Hamiltonian is then con- 
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structed, without any knowledge of the solution, to pro- 
duce a Rabi-like half-cycle which populates the marked 
state in a time scaling as N 1 / 2 , thereby exhibiting the 
quadratic speedup. An experimental realization of this 
analog Grover algorithm has been performed by NMR 
[111 ] in a setting where a quadrupolar coupling makes a 
spin 3/2 nucleus a two-qubit system (Af=4). 

Adiabatic processes offer many advantages, in partic- 
ular, the high degree of population transfer and some ro- 
bustness with respect to fluctuations of the control fields 
or imperfect knowledge of the model. Adiabatic versions 
of the time continuous Grover al gori thm have been con- 
structed in abstract form [l2|, fl3l Il4| . An ad hoc Hamil- 
tonian connects adiabatically the initial ground super- 
position to the marked state through an effective two- 
state avoided crossing. It has been shown [l4| that the 
transfer to the marked state in a time growing as N 1 / 2 
requires a specific time-dependent sweeping of the pa- 
rameter controlling the Hamiltonian. We have recently 
proposed a physical implementation of such an adiabatic 
search using three-state atoms trapped in a QED cavity 
and driven by laser fields [l5|. This scheme leads to an 
effective three-state dynamics which is closely related to 
stimulated Raman adiabatic passage (STIRAP) [l6j]. We 
have determined analytically the shapes of the adiabatic 
pulses that are required to lead to a quadratic speedup 
of the search. 

The present work contains a detailed description of 
the techniques and results announced in Ref. [15|]. We 
consider the more general case of search problems which 
have more than one solution. The dynamics of the col- 
lective marked, unmarked and excited states introduced 
below is solved exactly under local adiabatic conditions 
for which the scaling is determined analytically. Fur- 
thermore, we study numerically the robustness of the 
Grover search in this system using pulses which are eas- 
ily generated in practice (typically Gaussian pulses with 
plateaus). We identify the pertaining conditions which 
give rise to a speedup of the search which is only approx- 
imately quadratic in this case. 
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FIG. 1: Linkage pattern for the individual atoms. The un- 
marked atoms have two degenerate ground states \g) and \g'). 
The marked atoms feature the state \g') shifted. The laser of 
Rabi frequency Q' (resp. Q) is resonant with the g — e tran- 
sition for the marked (resp. unmarked) atoms. The cavity of 
Rabi frequency G is resonant with the g — e transition. 



The paper is organized as follows. In Sec. I, we de- 
scribe the cavity-laser-atom system, introduce the rele- 
vant collective states and derive the effective Hamilto- 
nian. Section III is devoted to the analysis of the adia- 
batic conditions leading to a Grover search with an ex- 
act or approximate quadratic speedup. The two types of 
pulses schemes are illustrated and discussed in Sec. IV 
while the conclusions are given in Sec. V. 



II. CAVITY-LASER-ATOM SYSTEM 

A. Description 

We use an ensemble of N identical three-level atoms 
trapped in a single-mode cavity of coupling frequency G, 
and driven by two lasers of Rabi frequencies and W . 
The key of the search process is that all the atoms have to 
be entangled initially as defined below. The atoms have 
a A configuration with two ground Zeeman states \g) and 
\g'), coupled to the excited state |e) by respectively the 
cavity and the two lasers (see Fig. 1). The states \g') 
of all the atoms are considered as the database. The 
state \g) stands here to couple all the atoms together by 
exchanging a single photon in the cavity. It allows us 
to consider the unmarked collective state, defined as the 
normalized sum of all the possible states \g'j) for all the 
unmarked atoms in state \g) except the j th in \g'). The 
state \g') of any marked atom is shifted, for instance by 
a magnetic field, such that the transition g' — e is res- 
onant with the laser Q' (resp. f2) for a marked atom 
(resp. unmarked atom). In this scheme, the magnetic 
field should be atomic-selective and the atoms should 
thus be fixed and form a register (of one or two dimen- 
sion) in the cavity. We start from the initial entangled 
state \g',0) = (1/y/W) $Zj =1 \g'j,ty featuring a collective 
superposition of both types of ground states. Such a 
state can be prepared for instance before the marking of 
the atom using the stimulated Raman adiabatic passage 
(STIRAP) technique, exactly as shown by Fleischhauer 



and coworkers to store single-photon quantum states [17j | . 
In their scheme, a single-photon wave packet enters the 
cavity through a mirror while the resonant laser Q is 
on and strong as !1 > Gy/~N, and the atoms are all in 
their ground state g). This generates the state denoted 
\g, 1). The laser is next switched off adiabatically while 
the photon is in the cavity, such that the population is 
transferred by STIRAP from the state \g, 1) to the collec- 
tive state ls',0). In this final state the photon is stored, 
shared among all the atoms. 

A related adiabatic process is constructed here in or- 
der to transfer adiabatically the initial entangled state 
\g',0) to a superposition of the marked states using an 
inverse fractional stimulated Raman adiabatic passage 
(if-STIRAP). This if-STIRAP is the time inversion of 
the fractional STIRAP (f-STIRAP) which transfer pop- 
ulation from a single state to a superposition of states 
[l8l |. The if-STIRAP allows one here to transfer the pop- 
ulation from the superposition \g\ 0) to the marked state. 



Model Hamiltonian 



The Hamiltonian describing the system of N atoms is 



M 



N 



(1) 



We consider a cavity mode of frequency to and coupling 
strength G together with two lasers of frequencies w, to— S 
and pulse shapes fi(t), f2'(t) which do not grow with N. 
The resonant driving provided by the cavity-laser-atom 
system is described by 

N 



V = uia^a + Ga Y"' \ej) (gj 



N 



ne wt + fiV< w - tf >*J WjWi I + h - c - ( 2 ) 

The full Hamiltonian H = Hq + V is block diagonal. Each 
block may be labeled by the number k of photons which 
are present in the cavity when the N atoms are in their 
ground state \g). The corresponding multipartite state 
\gi • • ■ 3jv) ® \k) is denoted \g, k). All the states which are 
coupled to \g, k) span a subspace whose projection oper- 
ator is Vk- The Hamiltonian may thus be rewritten as 
H = J2 k =o 'PkHVk- We shall focus on the block P 1 HV 1 
associated with a single photon in the cavity. 

The multipartite state \g, 1) is coupled by the cavity to 
any state \gi ■ ■ ■ gj-iejgj+i ■ • • ) ® |0) = \ej, 0) for which 
the atom j is in its excited state while the other atoms 
remain in their ground state \g). Each state \ej, 0) is cou- 
pled by the laser to the state \gi ■ ■ ■ gj-ig'^gj+i ■ ■ ■ gjy) <S> 
|0) = \g'j,0) where the atom j is in the ground state \g') 
while the other atoms are in the ground state \g). This 
coupling scheme is illustrated in Fig. 2 a, and the cor- 
responding projection operator, spanning this subspace 
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FIG. 2: (a) Coupling scheme in the cavity. The cavity G, 
laser Q, laser Q' Rabi frequency are featured by respectively 
thick, thin, and dashed arrows, (b) Equivalent scheme where 
the states 1^,0} and |e i7 0}), i = 1,M (resp. i = M+ 1,N) 
of frame (a) form the collective marked ground (Ig'^jQ)) and 
excited (|e m ,0)) states [resp. the collective unmarked ground 
(\g'u, 0}) and excited (|e u , 0}) states]. The effective cavity Rabi 
frequency to the collective marked (resp. unmarked) excited 



state is VMG (resp. y/N - MG). 
closer under H, reads 

N 

^ = E(^°>+^-o>)+?w 

3=1 



(3) 



C. Collective marked and unmarked dressed states 

In order to remove the oscillatory time dependence in- 
troduced by diagonal terms, we consider atomic states 
which are dressed by laser and cavity photons through 
the resonant transformation 

M N 

R = e- i5t ]T| 5 ;,o)( 5 ;,o|+ £ |^o><^,o| 

3=1 j=M+l 
N 

+ e-^ t \ 9 ,l){g,l\+e-^ t Y J \e j Me j A- (4) 

3=1 



This results in the new Hamiltonian 

H R = R\TiHVi)R - iR?R, 



(5) 



where the dot denotes time derivative. 

Among the marked atoms, none plays a privileged role. 
Similarly, the unmarked atoms are all equivalent. It is 
therefore appealing to treat them collectively. We con- 
sider the uniform superposition of marked (resp. un- 



M 



3=1 



/V 



K,0) = ^== £ W.0>, (6) 

v j=M+l 

and the uniform superposition of excited states associ- 
ated with the marked (resp. unmarked) atoms 



M 



,0) - ^Efeo). 



3=1 



Ie u ,0) = ^=L= £ |e,,0). (7) 

v j=M+l 

These states, together with \g,l), define a subspace which 
is relevant for the Grover search. Indeed, it can be shown 
that this subspace is closed under Hr, 

H r = VH r V + (I-V)H r (V 1 -V), (8) 

where 

V = V W m ,0) +V\BLfi) +^|e m ,0> +V\e n ,0) +V\g,l)- (9) 

We may thus restrict our attention to the Hamiltonian 

Hi = VHrV 

= G (VM\e m , 0) + VN- M\e u , 0)) (g, 1| 

+ S'| 5l ' n ,0)(e m ,0|+S| 3 ;,0)(e u ,0|+h.c, (10) 

where £ = fl + e~ m Vt' and £' = fi' + e M n. This Hamil- 
tonian is schematically represented in Fig. 2 b. 

D. Effective Hamiltonian 

The Hamiltonian H\ can be rewritten in a way which 
is prone to a further reduction. Indeed, notice that the 
first term of (|10|) actually features a linear combination 
of |e m , 0) and |e u , 0) which is nothing but the uniform su- 
perposition over all atoms ^7= X^=i I , 0) = |e, 0). This 
is expected since the state \g, 1) is coupled indifferently 
to any of the state | , 0} , marked or not, as is seen for 
instance in Fig. 2 a. By contrast, owing to the detuning, 
the g' — e transitions are not equivalent in the marked 
and unmarked cases. Hence, the last two terms of (fT0|) 
involve either |e m ,0) or |e u ,0). We may express Hi in 
terms of |e,0) and a state |ej_,0) which is orthogonal to 
the uniform superposition, 



|e,0) = y/f\e m ,0) + y/T=f\e u ,0), f 
|e±,0) = x/l - 7|e m ,0)-v7|eu,0). 



M 
~N 



(11) 
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FIG. 3: Coupling scheme in the A system corresponding to 
the effective Hamiltonian (|f 5|l with a fraction /(= M/N) of 
solution to the search problem. 



The result reads 



Hi — H e± + H e + H, 



e_L e ~\~ H-ee± j 



(12) 



where H e and H e± are coupled to each other only through 



H e±e and H ee± = Hl ±e , 



H e± =y/T=fTf |<4,0)(e ± ,0|-V7s K,0)(ex,0|+h.c. 
H e = y/NG |e,0)( 5 ,l| + h.c. 

(13) 



H e±e = ^' |<4,0}(e,0| + v/WS |^,0)(e,0 



As just mentioned, it is the Hamiltonian H e± which 
discriminates the marked and unmarked states whereas 
H e only features uniform superpositions over all the 
atoms. The subspace of interest is thus spanned by the 
projector 



Ve, = V 



(14) 



With the help of a unitary transformation T — exp(W) 
we can perturbatively block diagonalize the Hamiltonian 
H±. In this relevant subspace, the first correction to H e± 
is obtained by the partitioning technique, which is re- 
called in appendix, as \V e± [H ee± , eW] V e± ■ This quan- 
tity, featuring denominators which are the difference of 
the eigenvalues of H e± and H e , decreases at least as 
1/NG 2 (with additional decreasing contributions due to 
/). As a consequence, we shall neglect this correction 
together with the nonresonant components in £ and E'. 
We therefore consider the effective Hamiltonian 



H cS = y/T=frt |^,0)(ex,0|-v7n K,0)(e ± ,0|+h.c, 
which is schematically depicted in Fig. 



III. GROVER SEARCH 

A. Starting and final points 

Our aim is to transfer adiabatically the population 
from the initial state \g',0), 



1 N 



(16) 



which gives no privileged role to any of the N states | g'j , 0) 
to a final state which coincides with the collective marked 
state lg'^,0) in a time which scales as \J\jj where / = 
M/N is the fraction of solutions. The population transfer 
mechanism is most easily revealed in the basis of the 
instantaneous eigenstates of H c g (t) 



|0>(t) = cos d(t) \g , m ,0)~sm9(t) \g' n ,Q) 



(17) 



|±A)(t) = -^(sin0(i) \g' m> 0)+ cos 9(t) \g' u ,0) ± |ex,0)) 



pertaining to the eigenvalues and ±A(t) where 



A(t) 



mi 2 {t)+fw{t). 



(18) 



Note that |0) has no component on the collective excited 
states |ej_, 0) and is therefore a so-called dark state which 
is immune to loss by spontaneous emission (in contrast to 
the states | ± A)). The instantaneous angle 8(t) is defined 
through the relation 



tan0(t) 



'l-fil'jt) 



(19) 



Requiring the instantaneous eigenstate (| 1T[) to coincide 
at the initial time U with the uniform superposition 



|</,o> = V7ls' m ,o) + Vwi^,o) 



(20) 



and at the final time with the collective marked state 
| g' m , 0) entails that 



tan0(ii) 



1-/ 



/ 



tan 0(tf 



0. 



(21) 



This implies that the two pulses must be switched on 
simultaneously, fi'(ti) = £l(ti), and that the pulse f2' is 
to be turned off before f2. This process can be named 
inverse fractional STIRAP (if-STIRAP) since it allows 
the transfer by a STIRAP type process from a superpo- 
sition of states to a single state [lfjj]. In the adiabatic 
representation (JT7J) , the effective Hamiltonian reads 



Hg = A(| + A)(+A| + l-AX-AI) 
i€ 



V2 



(l + A)(0| 



A>(0| - h.c. 



(22) 



5 



4z tan 9. In the adiabatic regime, the 



where 9 = rrr a a 

l+tan^ 6 at 

transitions between instantaneous eigenstates are negli- 
gible. This will be achieved if the Hamiltonian varies 
sufficiently slowly in time so as to keep 9 -C A. On the 
other hand, we wish to control the process duration and, 
in particular, to prevent it from becoming arbitrary large. 
This can be achieved in several ways, two of which shall 
be considered explicitly. 



see from (|22|) that A appears as factor in the Hamiltonian 
Hj$. Hence, we can define a new time r(t) = J t dsA(s) 
so that the Hamiltonian H^/A is time independent. It 
follows that the survival probability amplitude of the 
state |0(t)) is 

,„/ m ( ^H^\.. .. l + £ 2 cosAr(i) . . 
(0(t)| exp -ir(t)-f- \0(U)) = —s U , (28) 



A 



1 



B. Local adiabatic conditions 

In order to control the process duration, as proposed 
in [l4| , w e choose to require 6 and A to be in a constant 
(small) ratio e at all times, independently of A^: 



= eA. 



(23) 



Given a laser pulse Q, this equation will allow us to de- 
termine the pulse Q' which is needed to remain in the 
instantaneous eigenstate |0)(t) with a significant prob- 
ability throughout the process, starting from the uni- 
form superposition |0)(ti) = \g' , 0), and ending up in the 
marked state |0)(tf) = \g' m , 0) after some time if — tj which 
achieves the optimal scaling with N. 

Indeed, let us first rewrite A defined in (|18|) in terms 
of tan 9 through (Til)]) as 



A = + tan 2 9Q. 



(24) 



We then obtain from (|23|) a differential equation for tan 9, 
i. e., for the ratio Q'/Q 



dtan# 



1 + tan 2 



.3/2 



(25) 



Its solution satisfying the initial condition (|2TJ) reads 



Ct(t) 



(26) 



sA(t)} 



where we define A(t) = f duil(u). Upon specifying that 
at time t{ the ratio of the pulses vanishes, one deduces 
from (J2BJ that sA(t f ) = - /)//. Expressing the 
total area of the pulse in terms of its peak amplitude 
f2o an d the process duration T, A(t{) — QqT , we finally 
obtain 



•l-f 



f 



(27) 



This expression shows that the search duration scales as 
J -1 / 2 for a peak amplitude which is independent of / = 
M/N. Equivalently, we can increase the peak amplitude 
as / _1 / 2 for a constant duration T. 

We now determine the population which can be 
reached by adiabatic passage. With the choice 



where A = y/1 + e 2 . The system therefore stays in the 
state |0(i)) with a probability satisfying at all times 



Po(t) > 



1 



(29) 



The non-adiabatic losses (to the states | ± A) which con- 
tain some component on the excited state |ej_)) are there- 
fore never larger than 4e 2 . Furthermore, since the states 
|0) and \g' m , 0) coincide at the final time, (f2"9")l also implies 
that the collective marked state ends up with a popula- 
tion of the order of 1 — 4e 2 . 

Returning to the diabatic representation, in which the 
state vector is denoted \4>(t)) , one can determine the full 
population dynamics of the collective states featured in 
the effective Hamiltonian (fT5|) . The population of the 
collective marked state is found to be 



P m (i) = \{g' m ,0\m\ 



T H.n Xr(t) sm9(t) + 1 + ^^W) C os6»(i)| 2 , 

(30) 



where 9(t) is given by p^|) and (PB]) . Recalling from (J^I 
that, at the final time, sin#(tf) = 0, one deduces that 



l + e 2 cosAT(t f )\\ /l-e 2x2 



1 + e 2 



> 



1 + e 2 



(31) 



The population P u (t) = \ (g' u , 0\<p(t)) | 2 of the collective 
unmarked state is given by an expression similar to (|30p 
where s'm9(t) and cos 9(t) are interchanged. In partic- 
ular, this entails that, at the final time, it is of order 

r-2 



P u (tf ) = t—^ sin2 Ar (*f 
1 + e z 



(32) 



Finally, the population of the collective excited state |ej_) 
is 



PeAt)= (1 f £2)2 (l-cosAr(t)) 2 , (33) 

which puts a bound on the decoherence present in the 
system since the population of this excited state is never 
larger than e 2 . 
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C. Tailored Gaussian pulses 

The local adiabatic conditions considered in the pre- 
ceding section allowed us to determine explicitly the scal- 
ing of the search duration and the population of the col- 
lective marked state. This is only a particular choice. 
The proposed scheme is robust in this respect. Indeed, 
as long as the initial and final conditions (|2T|) are satis- 
fied and one goes adiabatically from one to the other, the 
desired transfer of population can be achieved in a time 
which scales appropriately. In practical applications, it 
may be easier to consider Gaussian pulses possibly with 
plateaus. 

In order to satisfy the initial condition, the two pulses 
SI and SI' that are used in the current scheme are to 
be turned on simultaneously. The pulse SI' is a stan- 
dard Gaussian pulse Q' — Sloe"*-*/ 1 ^ . As for SI, it has 
a Gaussian switching on, and when the peak is reached, 
the pulse is kept at this value for a time aT before being 
switched off according to a Gaussian, 



15 



fi e- ( ' /T)2 

s> = { n a 



t > o 

< t < aT 
t > aT 



(34) 



We show numerically in the next section that with such 
pulses the search duration has a scaling which is close to 
optimal for a > 1.5. 
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FIG. 4: Analytical dynamics for the effective Hamiltonian 
(fT5|) with / = M/N — 3/8, e = 0.05 and Rabi frequencies 



given by a Gaussian profile f2(i) = Sloe 
(f26l) for fi'(i). Here Q T = ^/(T^ 
satisfy the conditions (|21|l and (|23|l 



-{t/ry 



and expression 



f)/f/ £ V^ m order to 
The final population 
transfer to the collective marked state is larger than 0.99. 
Top: Rabi frequencies. Bottom: Populations of the collective 
marked (Pm), unmarked (P u ) and excited {P e ±) states as a 
function of dimensionless time. 



IV. DISCUSSION 



The current adiabatic search scheme is robust in the 
sense that it leaves some flexibility on the pulses that are 
used. We presented explicitly two choices of pulses which 
will now be illustrated and compared. 



We first consider the case of (f2"3"|) which amounts to 
requiring local adiabatic conditions. The pulses are de- 
termined by (|26j) and displayed in the upper panel of Fig. 
II for / = M/N = 3/8, e = 0.05, and a Gaussian pulse SI 
of characteristic duration T. This value of e implies that 
the lower bound (|29p of the probability Pq to remain in 
the instantaneous eigenstate is 0.99. The population dy- 
namics of the collective marked, unmarked and excited 
states is given by the analytical expressions ([30|) - ([33|) and 
is depicted in the lower panel. As predicted by (|3Tj) . the 
transfer to the marked state is very efficient with a low 
transient population in the excited states stemming from 
the fact that the dynamics remains in the instantaneous 
decoherence-free eigenstate |0)(t) in the adiabatic limit. 
We remark that the choice (|23| , which leads to the seem- 
ingly complicated pulse relation (|2"6) . gives in practice a 
simple smooth bell-shaped pulse (see upper panel of Fig. 
HJ. Increasing / (i. e. increasing M for a given N) allows 
the reduction of the pulse area with the same efficiency, 
as shown by P?]). 



We now turn to the case of tailored Gaussian pulses. 
Figures [5] and [S] show the dynamics for such pulses with 
/ = M/N = 1/8 and / = 3/8 respectively. We have 
here chosen a plateau of duration 1.5T (a — 1.5). The 
pulse area of SI' has been chosen such that the popula- 
tion transfer to the marked state is approximately 0.99, 
i.e. as in the conditions of Fig. Figures [5] and [S] 
show similar features as [4j low transient population in 
the excited states, comparable pulse area for the same 
efficiency, reduction of the pulse area with the same effi- 
ciency for increasing M. We remark however two differ- 
ent properties: (i) For a similar efficiency, there is a larger 
transient population in the upper state for the cases of 
Figs. [5] and [SJ (ii) The population dynamics of the lat- 
ter cases do not exhibit the oscillations noticed in Fig. 
|4j Both features can be interpreted using superadiabatic 
basis that are better adapted to describe the dynamics 
fl9L |20T ] . This superadiabatic transport allows one to ex- 
plain the final revival of the transient "lost" population 
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FIG. 5: Numerical dynamics for the effective Hamiltonian 
(|15|) with / = M/N = 1/8 and Rabi frequencies chosen as 

Q'(t) = Sloe _( ' t '' T - ) and a Gaussian profile with a plateau for 
S}(t) as in (|34|) in order to satisfy (|21[) . Here QoT is such 
that the final population transfer to the marked state is 0.99. 
Top: Rabi frequencies. Bottom: Populations of the collective 
marked (Pm), unmarked (P u ) and excited (P e ±) states as a 
function of dimensionless time. 

in the upper state shown in Figs. [5] and [6] [21, 22]. Fur- 
thermore, the use of analytic pulses in the cases of Figs. 
[5] and [6] explain the non-oscillatory transfer. For the pre- 
ceding case, the local adiabatic control (|23|) induces the 
condition (|26[) which prevents in general the analyticity 
of one of the pulse (here W if is chosen analytic). In 
this non-analytic case, through non-adiabatic transitions 
in a superadiabatic basis (of order related to the order of 
the discontinuous derivative), population transfer occurs 
in an oscillatory manner as shown in Fig. [T] 

Figure [JJ displays the scaling of the search duration T 
for the scheme which uses the pulses of the the form (|34|) 
as shown in the upper panel of Fig. [5J The duration 
of the plateau has been kept to 1.5T (a — 1.5) and the 
quantity QqT has been determined numerically from the 
requirement that the population transfer to the collec- 
tive marked state reaches 0.99. We obtain a duration 
search scaling essentially as /~ ' 53 , which demonstrates 
the approximate quadratic speedup of the search. 

Finally, we investigate the role of the plateau in the 
efficiency of the search. The scaling exponent (3 of the 
search duration for various plateau durations aT is an- 
alyzed in Fig. [5J For each value of a, the power (3 is 
determined numerically by a linear fit of the log-log rep- 
resentation of fioT as a function of 1// (obtained in the 
same conditions as for Fig. [7]). Figure [S] shows that the 
power f3 decreases to 0.53 as the dynamics becomes more 
adiabatic until a w 1.5 before increasing slightly and sat- 



FIG. 6: Same as Fig. but with / = M/N = 3/8. 

4n 1 1 1 / 



3.6 



x 3.2 

a 

2.8- 



2.4 / 

1.6 2.4 3.2 4 4.8 

ln(1/f) 

FIG. 7: Scaling of QqT, with T the process duration, as a 
function of the fraction / of solution for the search with tai- 
lored Gaussian pulses as in Fig. [5] The plateau duration is 
fixed to 1.5 T and the final population of the collective marked 
state is 0.99. The natural logarithm of QoT as a function of 
the natural logarithm of 1// obtained numerically (dotted 
line) is well fitted by the straight line of slope 0.53. 



mating to 0.6 for a > 2.5. 



Before concluding, let us make some remarks concern- 
ing the experimental implementation of this adiabatic 
Grover search. Firstly, the scheme described here re- 
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FIG. 8: Power /3 of the scaling f~ p of Q T as a function of 
the plateau durations aT for tailored Gaussian pulses. The 
scaling exponent (3 is obtained as in Fig. [7] (see text for de- 
tails). 



quires to trap atoms in a cavity, a task which can be 
achieved for instance using a standing-wave dipole-force 
trap [23j]. As a realistic A atomic scheme, we can consider 
the typical 2 3 Si — 2 3 Po transition in metastable helium 
which is of linewidth T ~ 10 7 s _1 and Rabi frequency 
fl ~ lO 8 ^ s" 1 (with the intensity / in W/cm 2 ). In or- 
der to neglect spontaneous emission, as assumed in @, 
we require the condition (O T) 2 > TT. It is fulfilled 
when fl n T ^> 1 and il 3> T, which are well satisfied 
in practice, e.g., for Iq ~ 10 4 W/cm 2 and T ~ 10 ns. 
Finally, the feasibility of fixing the ratio of two pulses 
(here essentially required at early times in the local adi- 
abatic conditions approach) has been shown in [2H us- 
ing acousto-optical modulation of a cw laser in such a 
nanosecond regime. This step is not required in the other 
approach using Gaussian pulses with plateaus. One can 
mark the atoms by an ac Stark shift using a magnetic 
field, for instance from a focused laser beam for spatial 
resolution. 



V. CONCLUSIONS 

We have investigated theoretically a physical imple- 
mentation of the adiabatic Grover search using a cavity- 
laser-atom system and robust processes related to STI- 
RAP. The calculation has been conducted with pulses 
based on the constraint (f2"3"| that has allowed us to prove 
the scaling analytically. We have checked the robustness 
of the J -1 / 2 scaling by numerical simulations using other 
less restrictive adiabatic pulse shapes satisfying (|2"T|) . The 
case of Gaussian pulses with plateaus, which are routinely 



used in laboratories, has been studied. In particular, a 
scaling close to optimal is obtained when, after a Gaus- 
sian switching on, one of the two pulses is kept constant 
for a time equal to (or larger than) 3/2 of its character- 
istic time. Finally, the experimental implementation of 
the presented scheme has been briefly discussed. 
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APPENDIX A: PARTITIONING 

Let us consider a system whose Hilbert space is par- 
titioned into two orthogonal Hilbert spaces by means of 
the projection operators Va and V B = 1 — Va- The 
pertaining Hamiltonian can be written as 

H = H A + H B +H A B + H ba =H a + H b + eV, (Al) 

where e is a formal small (or more precisely ordering) pa- 
rameter. The purpose is to achieve a perturbative block- 
diagonalization of H. The leading contribution from the 
off-diagonal blocks to the diagonal ones can be extracted 
by the partitioning technique (25|. There is a unitary 
transformation exp(eVF) such that 

e- eW He eW = H A + H B + ^ [eV, eW] + e 3 V , (A2) 

with the property that [eV, eW] = V A [eV,eW}r A + 
V B [eV,eW]V B . The leading correction to the diagonal 
blocks is thus of order 2 while the remainder (e 3 V ) is of 
order 3. It has been shown that eW — eW B A — tW BA 
with eW B A given by 

tffM .. £ M<MM , (A3) 

where ) denotes the eigenvector of Ha corresponding 
to the eigenvalue \£ (and similarly for H B ). Recall that 
e is only formally a small parameter (it can be set to 1). 
However, note from (|A3|) that eW B A and the high order 
terms in (|A2[) are smaller, the larger the differences in 
the eigenvalues of Ha and H B . 

We now apply these results to the Hamiltonian (I12p . 
taking Hi for H, V e± for Ta and V e for V B . The domi- 
nant contribution from the off-diagonal blocks H e±e and 
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H ee , to H e is obtained from (IA2 



Hf} = H e± + ±V e± [zV,eW]V e± 

= H e± + ~ (H e±e eW ee± + h.c.) , (A4) 



where eW ee± is given by (jA3|) . The eigenvalues of H e are 
±\/iV G and the corresponding eigenvectors read 

\±^NG) ee -j={\e,0)±\g,l)). (A5) 

In the resonant approximation (E = ft + e~ ll5t ft' « 51 and 
S' = ft' + e 4(5 *ft w ft'), the eigenvectors of i? ei are 

|0) = cos9 |^,0) -sin0 K,0) (A6) 

|±A) = (sine |^,0) + cos^ |^,0)±|ex,0)), 

with the instantaneous angle 9 defined through tan 8 = 

-V(W)/Wn and 



(A7) 



To determine the leading correction to H e± given in (|A4[) , 
we note that -ff 6i e = H e±e V e — H e±e |e,0) (e,0|. Hence, 



we simply have to compute (e, 0| eW ee± , which after some 
algebra, reads 



(e,0|eW eR = 



A ( v7«' sin e + x/T" 7 ?^ cos 6>) 

NG 2 - A 2 
A^/I - /(ft 2 - ft' 2 ) cos ( 
(NG 2 - A 2 )n 



(ex, 01 



(ex,0| 



The correction featured in (IA4I) is then 



AvT^J(fi 2 - ft' 2 ) cos 



(A^G 2 - A 2 )ft 

x(V7n' |,4,o> + yi - 7^ | ff ;,o)) ( e± ,o|.(A8) 

Recall that, by definition of the problem, the constraint 
imposed in this Grover search is that the pulse envelopes 
ft and ft' do not increase with N, i. e., are at most of 
order A^°. It then follows from (|A7|) that A is at most 
of order A^°. Note that the additional dependency on A^ 
introduced through / may only decrease A. As a conse- 
quence, the correction (|A8[) to H e± decreases at least as 
1/NG 2 and may be discarded, giving rise to the effective 
Hamiltonian H c ff = H e defined in flE 
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